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We present an overview of the Hartree-Fock-Bogoliubov (HFB) theory 
of nucleonic superfluidity for finite nuclei. After introducing basic con- 
cepts related to pairing correlations, we show how the correlated pairs are 
incorporated into the HFB wave function. Thereafter, we present deriva- 
tion and structure of the HFB equations within the superfluid nuclear 
density functional formalism and discuss several aspects of the theory, 
including the unitarity of the Bogoliubov transformation in truncated 
single-particle and quasiparticle spaces, form of the pairing functional, 
structure of the HFB continuum, regularization and renormalization of 
pairing fields, and treatment of pairing in systems with odd particle 
numbers. 



1. Introduction 

Nucleonic pairing is a ubiquitous phenomenon underlying many aspects of 
structure and dynamics of atomic nuclei and extended nuclear matterP^l 
The crucial role of nucleonic superfluidity lies in its emergent nature. In- 
deed, while the correlation energy due to pairing is a small correction to 
the nuclear binding energy, the superfluid wave function represents an en- 
tirely different phase described by new quasiparticle degrees of freedom. 
The road to this new phase is associated with a phase transition connected 
with a symmetry breaking, and this underpins the nonperturbative nature 
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of pairing. 

Many facets of nucleonic superfluidity - including those related to phe- 
nomenology and theory of pairing - are discussed in this volumeP Here, 
we outline several aspects of nucleonic superfluidity within the framework 
of the nuclear density functional theory (DFT) . The main building blocks 
of nuclear DFT are the effective mean fields, often represented by local 
nucleonic densities and currents. When compared to the electronic DFT 
for the superconducting state^H^ the unique features of the nuclear variant 
are (i) the presence of two kinds of fcrmions, protons and neutrons, (ii) the 
absence of external potential, and (iii) the need for symmetry restoration 
in a finite self-bound system. In the context of pairing, nuclear superfluid 
DFT is a natural extension of the traditional BCS theory for electrons^ and 
nucleonsP and a tool of choice for describing complex, open-shell nuclei. 

At the heart of nuclear DFT lies the energy density functional (EDF). 
The requirement that the total energy be minimal under a variation of 
the densities leads to the Hartree-Fock-Bogoliubov (HFB; or Bogoliubov- 
de Gennes) equations. The quasiparticle vacuum associated with the HFB 
solution is a highly correlated state that allows a simple interpretation of 
various phenomena it the language of pairing mean fields and associated 
order parameters. 

This paper is organized as follows. Section [2] describes the essentials of 
the general pair-condensate state. The HFB theory is outlined in Sec. [3] 
The Bogoliubov sea, related to the quasiparticle-quasihole symmetry of the 
HFB Hamiltonian is discussed in Sec. |4]and Sec. [5] is devoted to the form 
of the nuclear pairing EDF. The quasiparticle energy spectrum of HFB 
contains both discrete bound states and continuum unbound states. The 
properties of the associated quasiparticle continuum are reviewed in Sec. [7] 
Section |8] describes the extension of the HFB formalism to odd-particle 
systems and quasiparticle blocking. Finally, conclusions are contained in 
Seel 

2. Basics of pairing correlations 

In quantum mechanics of finite many-fermion systems, pairing correlations 
are best described in terms of number operators = a+a^, where /i rep- 
resents any suitable set of single-particle (s.p.) quantum numbers. Thus 
we may have, e.g., /i = ha for plane waves of spin-^ particles (electrons); 
/i = rtJT for spin-i and isospin-i nucleons localized in space at position r; 
and n = n, £, j, m for fermions moving in a spherical potential well. Since 
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N„, the number operators are projective; hence, one can - at least 



in principle - devise an experiment that would project any quantum many- 
fermion state \^f) into its component with exactly one fermion occupying 
state ix. As the rules of quantum mechanics stipulate, any such individual 
measurement can only give or 1 (these are the eigenvalues of N^), whereas 
performing such measurements many times, one could experimentally de- 
termine the occupation probabilities = ( ^ 1 7V M | ^ ) - 

Along such lines, we can devise an experiment that would determine the 
simultaneous presence of two fermions in different orthogonal s.p. states \i 
and v. Since the corresponding number operators and N v commute, 
one can legitimately ask quantum mechanical questions about one-particle 
occupation probabilities and well as about the two-particle occu- 

pation. The latter one reflects the simultaneous presence of two fermions 
in state |^}: v^ v = {^iN^N^ty). In this way, one can experimentally 
determine the pairing correlation between states /x and v as the excess 
probability 



of finding two fermions simultaneously over that of finding them in in- 
dependent, or sequential, measurements. Such a definition of pairing is 
independent of its coherence, collectivity, nature of quasiparticles, symme- 
try breaking, thermodynamic limit, or many other notions that are often 
associated with the phenomenon of pairing. In terms of occupations, pair- 
ing can be viewed as a measurable property of any quantum many-fermion 
state. 

Obviously, no pairing correlations are present in a quantum state that 
is an eigenstate of jV m or N„, such as the Slater determinant. The beauty of 
the BCS ansatz is in providing us with a model ./V-fermion state, in which 
pairing correlations are explicitly incorporated: 



where the summation fi > runs over the representatives of pairs (ft, fi) of 
s.p. states (that is, any one state of the pair is included in the sum, but not 
both), = z p are real positive numbers, = — s M are arbitrary complex 
phase factors, and M m is the overall normalization factor. 

It now becomes a matter of technical convenience to employ a particle- 
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number mixed state, 



l*>=^ E Mf ^ Ew^ lO), (3) 

JV=0,2,4,... V 1 '' \jU>0 / 

in which the pairing correlations ([T]) are: 

z 2 1 
= vlulSp, u for ^ = /* 2 and = — — . (4) 

I -\- Z^ 1 -f- Z v 

In terms of the s.p. occupations, the state |$) assumes the standard BCS 
form: 

i$>=n(^+w4<)i°>- (5) 

M>0 

In this many-fermion state, the s.p. states p and /x are paired, that is, 
|<fr) can be viewed as a pair-condensate. For the pairing correlation 

Pp,n © equals 1/4; in fact, in this state, it is twice more likely to find a 
pair of fermions (w|n = l/2) than to find these two fermions independently 
(u? 1^=1/4). For particle- number conserving states ([2]), the occupation 
numbers can be calculated numerically; qualitatively the results are fairly 
similar, especially for large numbers of particles. 

At this point, we note that the most general pair-condensate state (|3]) 
has the form of the Thouless state, 



|0>, (6) 

in which pairs (p, /i) do not appear explicitly. However, there always exists 
a unitary transformation Uo of the antisymmetric matrix Z that brings it 
to the canonical form (UqZ*Uq ) Vfl — s^z^d^ (the Bloch-Messiah-Zumino 
theorenP^ . Therefore, pairs are present in any arbitrary Thouless state 
(the so-called canonical pairs), and they can be made explicitly visible by 
a simple basis transformation. 

The canonical pairs exist independently of any symmetry of the Thou- 
less state. In the particular case of a time-reversal-symmetric state, 
T|$) = |$), they can be associated with the time-reversed s.p. states, 
jl = p,. The ground-states of even-even nuclei can be described in this man- 
ner. However, the appearance of pairing phase does not hinge on this partic- 
ular symmetry - states in rotating nuclei (in which time-reversal symmetry 
is manifestly broken) can also be paired. In this latter case the canonical 
states are less useful, because they cannot be directly associated with the 
eigenstates of the HFB Hamiltonian. 
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3. Hartree-Fock-Bogoliubov theory 

The simplest route to the HFB theory is to employ the variational prin- 
ciple to a two-body Hamiltonian using Thouless states © as trial wave 
functions. The variation of the average energy with respect to the antisym- 
metric matrix Z results in the HFB equation in the matrix representation, 
WA = UE, or explicitly, 

fT + Y A UP V*\_ fU V*\ (E \ 

^ -A* -t* - r* ) \ v u* ) ~ \ v u* ) \o -e ) ' i, ] 

where T M „ is the matrix of the one-body kinetic energy, T ^ — 

Y*n> V ' Vw';w'Pv'im' and A^' = ^E W 'V;»'V are tne so-called 
particle-hole and particle-particle mean fields, respectively, obtained by 
averaging two-body matrix elements V^-^v' with respect to the density 
matrix p v i^ — ($|a+,a„' |$) and pairing tensor n vv i = {§\a v i a v \$>) , and E 
is the diagonal matrix of quasiparticle energies.^ 

The matrices TL and U are referred to as the HFB Hamiltonian and 
Bogoliubov transformation, respectively, and columns of U (eigenstates of 
H) are vectors of quasiparticle states. The HFB equation (J7J possesses the 
quasiparticle-quasihole symmetry. Namely, for each quasiparticle state Xa 
(the a-th column of U) and energy E a there exists a quasihole state 4> a of 
opposite energy — E a , 





4>a = 








\u; a ) 



(8) 



That is, the spectrum of % is composed of pairs of states with opposite 
energies. In most cases, the lowest total energy is obtained by using the 
eigenstates with E a > as quasiparticles \a and those with E a < as 
quasiholes <j> a , that is, by occupying the negative-energy eigenstates. States 
Xa and 4> a can usually be related through a self-consistent discrete symme- 
try, such as time reversal, signature, or simplex E2HE] 

The HFB equation (JTJ) is also valid in a more general case, when the 
total energy is not equal to the average of any many-body Hamiltonian. 
Within the DFT, it stems from the minimization of the binding energy 
given by an EDF £(p, k, k*), subject to the condition of the generalized 
density matrix being projective, that is, TZ 2 = 1Z for 

/ p k \ fV*V T V*U T \ fV*\( V T U T ) V- - - + 

(9) 
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In this case, the mean fields are obtained as functional derivatives of EDF: 
= d£jdp vil and = dS/dn* ,. As is the case in DFT, densities 

(here the density matrix and pairing tensor) become the fundamental de- 
grees of freedom, whereas the state |$) acquires the meaning of an auxiliary 
entity (the Kohn-Sham state^). Indeed, for any arbitrary generalized den- 
sity matrix 1Z one can always find the corresponding state |<I>). For 
that, one determines the Bogoliubov transformation U as the matrix of its 

eigenvectors, VIA = U ( jj the Thouless state |$) ([6]) corresponds to 

Z = VU^ 1 . Consequently, the paired state |$) of DFT is not interpreted 
as a wave function of the system - it only serves as a model for determining 
one-body densities. Nonetheless, these densities are interpreted as those 
associated with the (unknown) exact eigenstate of the system. 

Unrestricted variations of the EDF are not meaningful. Indeed, since 
Thouless states © are mixtures of components with different particle num- 
bers, absolute minima will usually correspond to average particle numbers 
that are unrelated to those one would like to describe. In particular, for self- 
bound systems governed by attractive two-body forces (nuclei), by adding 
more and more particles one could infinitely decrease the total energy of 
the system. Therefore, only constrained variations make sense, that is, one 
has to minimize not the total energy £ (p, k, k*), but the so-called Routhian, 
£ '(p, k, k*) — £(p, k, k*) + C(p), where C is a suitably chosen penalty func- 
tional, ensuring that the minimum appears at prescribed average values of 
one-body operators. In particular, the average total number of particles 
can be constrained by C(p) = — A($|A r |$) = — ATr(p) (linear constraint) or 
C(p) = Cjsr[Ir(p) — iV ] 2 (quadratic constraint) 11SU2I where A becomes the 
Fermi energy corresponding to No fermions. 

For different systems and for different applications, various constraints 
C{p) can be implemented; for example, in nuclei one can simultaneously 
constrain numbers of protons and neutrons, as well as multipole moments 
of matter or charge distributions. When the total energy is a concave func- 
tion of relevant one-body average values, quadratic constraints are manda- 
£ 0r yji6ii7] rp-^g minimization of £'(p, k, k*) requires solving the HFB equa- 
tion for the quasiparticle Routhian % , which, for the simplest case of the 

constraint on the total particle number, reads H' = % — A f J ^ 

Finally, let us mention that in the coordinate space-spin(-isospin) rep- 
resentation, the HFB equation acquires particularly interesting form, 
which in condensed matter and atomic literature is called Bogoliubov-de 
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Gennes equation^ In the coordinate representation, quasiparticle vectors 
become two-component wave functions, which - in finite systems - acquire 
specific asymptotic properties'^ determining the asymptotic behavior of 
local densities. The quasiparticle energy spectrum of HFB contains dis- 
crete bound states, resonances, and non-resonant continuum states. As 
illustrated in Fig. [IJ the bound HFB solutions exist only in the energy re- 
gion \Ei\ —A. The quasiparticle continuum with \Ei\ > —A consists of 
non-resonant continuum and quasiparticle resonances, see Sec. [6] 
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Fig. 1. Quasiparticle spectrum of the HFB Hamiltonian. The bound states exist in 
the energy region ^ —A, where A is the chemical potential (always negative for a 
particle-bound system). 



4. Beneath the bottom of the Bogoliubov sea 

Because of the quasiparticle-quasihole symmetry (Sec. [3]), the spectrum of 
the HFB Hamiltonian contains as many negative as positive eigenvalues. 
Therefore, the HFB equation ([7]) constitutes an eigenvalue problem for the 
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operator unbounded from below, and the HFB spectrum extends from mi- 
nus to plus infinity, see Fig. [TJ Moreover, one-body densities (density ma- 
trix and pairing tensor) are given by infinite sums over all negative-energy 
(quasihole) states, cf. Eq. ((9J. In analogy to the Fermi sea of occupied 
states, which appears in the Hartree-Fock (HF) theory, we call the set of 
quasihole states the Bogoliubov sea. Note again that the Fermi sea extends 
over a finite range of energies - from the bottom of the HF potential up 
to the Fermi energy - whereas the Bogoliubov sea is infinitely deep, in a 
nice analogy with the sea of one-electron states of the relativistic Dirac 
equation. 

In practice, since infinite sums over the Bogoliubov sea cannot be carried 
out, the number of pairing-active states must be truncated. Two different 
ways of achieving this goal are most often implemented, namely, solution 
of the HFB equations in a finite s.p. space (e.g., the so-called two-basis 
method^) and truncation of the summation in the quasiparticle space. 
The second method correspond to creating an artificial finite bottom of the 
Bogoliubov see. In this section we discuss the consequences of neglecting 
the quasihole states that are below the bottom of that sea! 23 * 24 ! 

The main problem concerns the calculation of the pairing tensor, which 
is the sum of products of upper and lower components of quasihole states, 
cf. Eq. ((9]). When this sum is performed over the infinite complete set of 
quasiparticle states, the resulting pairing tensor is antisymmetric, while for 
truncated sums it may acquire a symmetric part. Usually the symmetric 
component is small}^ hence, can be neglected. However, its very existence 
means that the many-fermion state which would have had such a pairing 
tensor simply does not exist. The smallness of the symmetric part can be 
deceiving, because the symmetric pairing tensor corresponds to a many- 
boson system. Consequently, appearance of the symmetric component im- 
plies the violation of the Pauli principle. This is a potentially dangerous 
situation - within a variational theory one should avoid the boson sector 
whose ground state is way below the fermionic ground state. 

A solution to this problem,^ discussed below, consists in marrying the 
two truncation methods mentioned above. That is, we shall first use the 
quasiparticle truncation method to define the appropriate s.p. cutoff, and 
then the HFB equations are solved in this truncated space, leading to a 
perfectly antisymmetric pairing tensor. 

Let us consider the case of truncated summations over the Bogoliubov 
sea and assume that we have kept only K quasihole states. In order to main- 
tain the quasiparticle-quasihole symmetry, we apply the same truncation 



March 8, 2013 0:51 World Scientific Review Volume - 9in x 6in BCS50-28 



Hartree-Fock-Bogoliubov solution of the pairing Hamiltonian in finite nuclei 9 

to quasiparticle and quasihole states, that is, we also keep K quasiparticlc 
partner states. This is convenient, and always possible, because the quasi- 
particle (unoccupied) states do not impact HFB densities. Then, matrix U 
becomes rectangular - it has less columns (2K) than rows. Since all kept 
quasiparticle and quasihole states are orthonormal, we still have U + U = 1. 
However, since now the quasiparticle space is not complete, U is not any- 
more unitary, and the product UU + = V is not equal to unity. All what 
remains is the hermitian and projective property of V: 

V + = P, V 2 = V. (10) 

Since TrT 5 = 2K, V has exactly 2K eigenvalues equal to 1. 
In its explicit form, matrix V reads: 

_ / p Q \ _ ( uu+ + v*v T , uv+ + v*u T \ 

~ \Q* P* ) ~ \VU+ + U*V T , VV+ + U*U T ) ■ 
In terms of P and Q, Eqs. (fTU)) can be written as: 

P+ = P, Q T = Q and P 2 + QQ+ = P, PQ + QP* = Q, 
where TrP = K. 

Properties of matrices P and Q can be most easily discussed in the 
particle basis that diagonalizes P. Suppose that column / is an eigenvector 
of P with eigenvalue p, that is, Pf — pf. From Eqs. {H} it follows that p 
must be between and 1. Moreover, if Qf* is not equal to zero, then Qf* 
is an eigenvector of P with eigenvalue 1—p, that is, P(Q/*) = (1 —p)Qf*. 
Conversely, if Qf* is equal to zero, than p=0 or p=l. 

Altogether, the spectrum of P can be divided into three regions: (i) 
i states with p v =l, where all matrix elements of Q vanish, Q UU '—0, (ii) 
2k states with 0<p„<l, where eigenvectors are arranged in pairs pc,=l— p v 
such that the only non- vanishing matrix elements of Q are 

Qua = Quu = Qu= VPJl ~Pv), (13) 

and (iii) states with p^—0, where again all matrix elements of Q vanish. In 
practical calculations of solving the HFB equation in infinite-dimensional 
quasiparticle spaces (like the coordinate representation), the first region 
almost never appears (i — 0), and then k — K. However, for truncated 
quasiparticle spaces, the third region always exists and contains the null 
space of P. 

We now see that when K quasiparticles are included in the quasiparticle 
space and i = 0, in the particle space there appears a basis of 2K s.p. states, 



(11) 
(12) 
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which we call natural basis. Each state in the first half of the natural basis 
has its partner in the second half. By ordering the eigenvalues of P and 
neglecting the zero eigenvalues for v > 2K, we can write matrices P and Q 
in a general form: 



(pi . 

P 2- 











■Pk 
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I-PK 
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Q = 



.. 
.. 



... 







■ -qK 





q 2 



q K 



















... l-p 2 ... o<2 2 ...ooo...o 
\ 0...0 ...1-pJ \ ?1 0...0 0...0 / 

(14) 

If there appear i > states with p u = 1, the number of paired states 
decreases to 2k = 2K — 2i and the size of the natural basis decreases to 
i + 2k = 2K - i. 

The HFB equations can now be solved in the finite natural basis, where- 
upon the pairing tensor becomes exactly antisymmetric and the dangerous 
violations of the Pauli principle are removed exactly.^ The advantage of 
this method is in the fact that the truncated s.p. space is not arbitrarily 
cut but it is adjusted to the truncated quasiparticle space. 



5. Pairing functional 

The form of the most general pairing EDF that is quadratic in local isoscalar 
and isovector densities has been discussed in Refs.^^Sl Because of the lack 
of nuclear observables that can constrain coupling constants of this general 
pairing functional, current realizations are much simpler. A commonly 
used effective pairing interaction is the zero-range pairing force with the 
density-dependent form factorPStHl 

Po(r)~ 



/pair(r) = V U + x P 



V- 



(l+x 3 P a ) 



(15) 



where P a is the usual spin-exchange operator and /9o=0.16fm -3 . When 
only the isovector pairing is studied, the exchange parameters xq and X3 
are usually set to 0. However, in the general case of coexisting isoscalar 
and isovector pairing correlations, nonzero values of xq and 23 have to be 
used. Pairing interactions corresponding to 77=0, 0.5, and 1, are usually 
referred to as volume, mixed, and surface pairing, respectivelyP^lES] The 
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volume pairing interaction acts primarily inside the nuclear volume while 
the surface pairing generates pairing fields peaked around or outside the 
nuclear surface. 

Another form of density dependence has been suggested in RefP^ and 
successfully applied^ to explain odd-even effects in charge radii. As dis- 
cussed in Refs.J^i^ different assumptions about the density dependence 
may result in differences of pairing fields in very neutron rich nuclei. How- 
ever, the results of the global survey^ suggest that - albeit there is a slight 
favoring of the surface interaction - one cannot reliably extract the den- 
sity dependence of the effective pairing interaction (j!5[) from the currently 
available experimental odd-even mass differences, limited to nuclei with a 
modest neutron excesses (see also RefsiMMM]), 

A timely question, related to the density dependence, is whether there 
is an effective isospin dependence of the pairing interaction. The global 
survejfUD f odd-even staggering of binding energy indicates that the effec- 
tive pairing strength Vq for protons is larger than for neutrons, and the 

recent large-scale optimizations of the nuclear EDF are consistent with this 
nn ding|M142] This 

can be attributed to the isospin-dependent contribution 
to pairing from the Coulomb interactiorP^SH or to induced pairing due to 
the coupling to collective excitations) 46 * 47 * To account for those effects, an 
extended density dependence has been proposecP^ES! that involves the local 
isovector density pi(r). 

Little is known about the isoscalar pairing functional. The local isovec- 
tor pairing potentia l 25 * 26 ! U(r) is proportional to the isovector pair density 
p whereas the isoscalar pairing potential 2o(t*) is a vector proportional to 
the isoscalar- vector pairing spin density Sq. Then, the isoscalar pairing 
field, 

ho(r) = S • o- oc s • <r (16) 

is the projection of the quasiparticle's spin on the proton-neutron pairing 
field. Physically, p represents the density of S=Q, neutron-neutron, proton- 
proton, and proton-neutron pairs, whereas the vector field Sq describes 
the spin distribution of 5=1 pn pairs (that is, it contains all magnetic 
components of S—l pn pairing field). 

Symmetries of the isoscalar pairing mean-fields have been studied in 
detail in RefP^ As an example, lines of the solenoidal field Sq - present in 
the generalized pairing theory that mixes proton and neutron orbits - are 
schematically shown in Fig. [2] It is interesting to note that for the geometry 
of Fig. [2j the third component s 0z , associated with the M=Q isoscalar 
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pairing field vanishes. That is, the solenoidal pairing field is created by the 
two components with M=±l. One can thus conclude that the assumption 
of axial symmetry, or signature, does not preclude the existence of isoscalar 
pairing. 
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Fig. 2. (a) Schematic illustration of the isoscalar vector pairing field so in the case of 
conserved axial and mirror symmetries!^ The field is solenoidal, with vanishing third 
component, (b) Under rotation around the third (symmetry) axis, the field at the point 
r& is transformed to the position rg. Likewise, under TZ X , the field is transformed to 
tq. While neither of these operations leave the individual vector so( r >t) invariant, the 
field as whole does not change. 



6. Pairing and the HFB continuum 

The structure of HFB continuum indicated schematically in Fig.[T]has been 
a subject of many works jH} l 21 l 51 l tHl Within the real-energy HFB frame- 
work, the proper theoretical treatment of the HFB continuum is fairly so- 
phisticated since the scattering boundary conditions must be met. One 
way of tackling this problem is the coordinate-space Green's function tech- 
niqueP3Ml56| jf ^ ne outgoing boundary conditions are imposed, the un- 
bound HFB eigenstates have complex energies; their imaginary parts are 
related to the particle widthPSl The complex-energy spherical HFB prob- 
lem has been formulated and implemented within the Gamow HFB (GHFB) 
approach of RefP^ 

In addition to the methods that directly employ proper asymptotic 
boundary conditions for unbound HFB eigenstates, the quasiparticle con- 
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tinuum of HFB can be approximately treated by means of a discretization 
method. The commonly used approach is to impose the box boundary 
conditions in the coordinate-space calculations l 19 l 21 l 57 l 59 < MI i n which wave 
functions are spanned by a basis of orthonormal functions defined on a 
lattice in coordinate space and enforced to be zero at box boundaries. In 
this way, referred to as the C 2 discretization, quasiparticle continuum of 
HFB is represented by a finite number of box states. It has been demon- 
strated by explicit calculations for weakly bound nuclei 5 - 2 ^- 57 * 58 ! that such 
a box discretization is accurate when compared to the exact results. Al- 
ternatively, diagonalizing the HFB matrix in the Poschl-Teller-Ginocchio 
basis^ or Woods-Saxon basis^flM] turned out to be an efficient way to ac- 
count for the continuum effects. Finally, quasiparticle continuum can be 
effectively discretized by solving the HFB problem by means of expansion 
in a harmonic oscillator (HO) or transformed HO (THO) basis P^^D As far 
as the description of nonlocalized HFB states is concerned, the coordinate- 
space method is superior over the HO expansion method, as the HO basis 
states are always localized. Consequently, the discretized representation of 
the quasiparticle continuum is different in coordinate-space and HO basis- 
expansion approaches.^ 

Among the quasiparticle resonances, the deep-hole states play a dis- 
tinct role. In the absence of pairing, a deep-hole excitation with energy 
Ei > corresponds to an occupied HF state with energy e$ = — Ei + A. 
If pairing is present, it generates a coupling of this state with unbound 
particle states with £j w Ei + A that gives rise to a quasiparticle reso- 
nance with a finite width.^ 21 * 72 ! Quasiparticle resonance widths can be 
directly calculated with a high precision using coordinate-space Green's 
function techniqu e 120155 1 56 1 and GHFBpSl For approaches based on the C? 
discretization, several approximate methods have been developed to deal 
with HFB resonances. The modified stabilization method based on box 
solution^ 57 1 73174 1 has been used to obtain precisely the resonance energy and 
widths. Based on the box solutions, the HFB resonances are expected to 
be localized solutions with energies weakly affected by changes of the box 
size. The stabilization method allows to obtain the resonance parameters 
from the box-size dependence of quasiparticle eigenvalues. 

Besides the stabilization method, a straightforward smoothing and fit- 
ting method that utilizes the density of box states has been successfully 
used. In this technique, resonance parameters are obtained by fitting 
the smoothed occupation numbers obtained from the dense spectrum dis- 
cretized HFB solutions. Figure [3] displays occupation probabilities vf for 
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Fig. 3. Occupation numbers of the discretized neutron quasiparticle continuum states 
calculated for 70 Zn in Ref.ES The corresponding Breit-Wigner envelopes are indicated 
by dashed lines. The — A n threshold is marked by a dotted line. 



the discretized neutron quasiparticle states in T Zn as a function of quasi- 
particle energy £j. To extract resonance parameters from the discrete dis- 
tribution of vf , one can first smooth it using a Lorentzian shape function 
and then perform a fit using a Breit-Wigner function.^ 

Various ways of computing widths of high-energy deep-hole states have 
been compared in Refp3 By comparing with the exact GHFB results, it 
has been concluded that the stabilization method works fairly well for all 
HFB resonances, except for the very narrow ones. The smoothing-fitting 
method is also very effective and can easily be extended to the deformed 
case. The perturbative Fermi golden rule^l has been found to be unreli- 
able for calculating widths of neutron resonances. (For more discussion on 
limitations of the perturbative treatment, see RefP^). 

Pairing correlations can profoundly modify properties of the system in 
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drip line nuclei due to the presence of the vast continuum space avail- 
able for pair scattering. One example is the appearance of the pairing- 
antihalo effect jU MEH EH i n which pairing correlations in the weakly-bound 
even-particle system change the asymptotic behavior of particle density 
thus reducing its radial extension. While neutron radii of even-even nuclei 
are expected to locally increase when approaching the two-neutron drip 
lin e! 37 * 65 * 72 * 77 * 79 ! the size of the resulting halo is fairly modest, especially 
when compared with spatial extensions of neighboring odd-neutron sys- 
tems. 

Pairing correlations impact the limits of the nuclear existence: the odd- 
even staggering of the nuclear binding energy does result in the shift be- 
tween one-neutron and two-neutron drip lines. The pairing coupling to 
the positive-energy states is an additional factor influencing the nuclear 
binding? 20 ! 21 ! In particular, because of strong coupling to the neutron con- 
tinuum, the neutron chemical potential may be significantly lowered thus 
extending the range of bound nuclei, and this effect is expected to depend 
on the character of pairing interaction. For more discussion on the impact 
of continuum on quasiparticle occupations, emergence of bound canonical 
HFB states from the continuum, and contributions of nonresonant contin- 
uum to the localized ground state in dripline nuclei, see Refs! 2 ^ 7 * 56 * 68 * 72 * ^ 



7. Regularization of the local pairing interaction 

As discussed in Sec. [5l in many HFB applications, pairing interaction is 
often assumed to be in the form of the zero-range, density-dependent force. 
Calculations using the contact interaction are numerically simpler, but the 
pairing gap diverges when the dimension of the pairing-active space in- 
creases for a fixed strength of the interaction. In roots of this problem is 
the ultraviolet divergence of abnormal density for zero-range pairing inter- 
action™™ 



h(v)M*{r) 

p(r - x/2, r + x/2) - 



(17) 



Consequently, in practical calculations, one has to apply a cutoff proce- 
dure to truncate the pairing-active space of s.p. states,^ 21 * 37 1 and the 
pairing strength has to be readjusted accordingly. Thus the energy cut- 
off and the pairing strength together define the pairing interaction, and 
this definition can be understood as a phenomenological introduction of 
finite ranged 82 ! 83 ! Such a sharp cut-off regularization is performed in the 
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spirit of the effective field theory, whereupon contact interactions are used 
to describe low-energy phenomena while the coupling constants are read- 
justed for any given energy cutoff to account for high energy effects. It has 
been shown that by an appropriate renormalization the pairing strength for 
each value of the cutoff energy, one practically eliminates the dependence 
of various observables on the cutoff parameter! 21 ! 54 ! Figure 2] illustrates the 
procedure for the total energy in the tin isotopes. While for a fixed pair- 
ing strength total energies depend significantly on the cut-off energy (top) , 
for a fixed pairing gap the changes obtained with renormalized interactions 
(bottom) are very small indeed. 



> 



< 



/ 



Sn 



v„ =-160 MeV 



AE mI ^40 MeV 1 



/ SkP 8 



V„ renormalized ~\Jl!l]^ 



AE m =40 



MeV | — ^-\^/?' 



50 70 90 110 130 

Neutron Number 

Fig. 4. Total energies in the tin isotopes calculated within the HFB+SkP 5 modelPS 
Top panel shows the results for the fixed interaction strength Vq and for several cut-off 
energies A_E max added to the usual £j-dependent cut-off energy E max P^ Bottom panel 
shows similar results when the values of Vb are renormalized to keep the average neutron 
pairing gap in 120 Sn the same for each A_E max . 

The cutoff energy dependence of the pairing strength can also be han- 
dled by means of a regularization scheme by defining the regularized local 
abnormal density! 54 ' 8 " ' 81 ' 84 ! ^ 



p r (v) = lim [p(r - x/2, r + x/2) - f(r, x)] 

x— >0 



(18) 



where / is a regularization countcrtcrm, which removes the divergence (If 7[) 
at x = 0. For cutoff energies high enough, one can express / through the s.p. 
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Green's function at the Fermi level, G(r + x/2, r — x/2), which also exhibits 
a 1/x divergence. In practical calculations, one cane use the Thomas- 
Fermi (TF) approximation for the local s.p. Greens function; this approach 
has been used with success for a description of spherical and deformed 
nuclei.^ 54 * 86 * 87 ! As demonstrated in RefP^ the differences between pairing 
renormalization and regularization procedures are rather small. 

A combination of the renormalization and regularization methods de- 
scribed above is the hybrid technique^ based on the TF approximation to 
the non-resonant HFB continuum! 88 ! 89 ! This approach is of great practi- 
cal interest as it makes it possible to carry out calculations in wide pair- 
ing windows and very large coordinate spaces. In the hybrid method, the 
high-energy continuum above the cutoff energy E c is divided into the non- 
resonant part and deep-hole states. While deep-hole states have to be 
treated separately, the non-resonant continuum contribution to HFB den- 
sities and fields can be integrated out by means of the TF approximation. 
The choice of the cutoff E c is determined by positions of deep-hole levels}^ 
this information can be obtained by solving the HF problem. 

8. Pairing in odd-mass nuclei 

The zero-quasiparticle HFB state representing the lowest configuration 
for a system with even number of fermions, corresponds to a filled sea of 
Bogoliubov quasiholes with negative quasiparticle energies, see Fig. [1] In 
a one-quasiparticle state representing a state in an odd nucleus, a positive- 
energy quasiparticle state a is occupied and its conjugated quasihole part- 
ner is empty. The corresponding wave function can be written as 



which depends on the quasiparticle state \a ® • Density matrix and pairing 
tensor of state (1191) can be obtained by exchanging in U columns correspond- 
ing to the quasiparticle and quasihole states, Xa and (j> a . The corresponding 
density matrix reads explicitly, 




(19) 



where a+ is the quasiparticle creation operator, 




(20) 




(21) 
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and similar holds for the pairing tensor. After the column replacement, 
matrix U^l of one-quasiparticle state becomes singular and has null space 
of dimensions one. Hence, the occupation number of one of the s.p. states 
equals to 1 . This fact is at the origin of the name "blocked states" attributed 
to one-quasiparticle states (fT9|) . These states contain fully occupied s.p. 
states that do not contribute to pairing field .USHUl 

The blocking can also be implemented, for some configurations, by in- 
troducing two chemical potentials for different supcrfluid components (two- 
Fermi level approach, 2FLA^^ As demonstrated in Ref., 9 ^ such proce- 
dure is equivalent to applying a one-body, time-odd field that changes the 
particle-number parity of the underlying quasiparticle vacuum. For po- 
larized Fermi systems, in which no additional degeneracy of quasiparticle 
levels is present beyond the Kramers degeneracy, the 2FLA is equivalent to 
one-dimensional, non-collective rotational cranking. 

When describing properties of odd-mass nuclei, one selects the lowest 
quasiparticle excitations E a and carries out the self-consistent procedure 
based on these blocked candidates (Tl9l) . Naturally, one must adopt a pre- 
scription to be able to determine, at each iteration, the index a of the 
quasiparticle state to be blocked.^ Such a unique identification can be 
done by means of, e.g., the overlap method of Ref! 9 ^ After the HFB itera- 
tions are converged for each blocked candidate, the state corresponding to 
the lowest energy is taken as the ground state of an odd-mass nucleus, and 
the remaining ones are approximations of the excited states. A similar pro- 
cedure can be applied to many-quasiparticle states, e.g., two-quasiparticle 
states in even-even and odd-odd nuclei, three-quasiparticle excited states 
in odd-mass nuclei, and so on. 

The state (fT9|) represents an odd-Fermi system that carries nonzero 
angular momentum; hence, it breaks the time reversal symmetry. If the 
time reversal symmetry is enforced, additional approximations have to be 
applied based on the Kramers degeneracy. One of them is the equal filling- 
approximation (efa), in which the degenerate time-reversed states Xa and 
Xa are assumed to enter the density matrix and pairing tensor with the 
same weights!- 94 * 99 ^ For instance, the blocked density matrix of efa reads: 

pi a J' EFA = (y*v T )^ + \ {u^u; a - v; a v va + u, a u: & - v; & v v& ) . (22) 

It has been showrPEO that the efa and the exact blocking are both strictly 
equivalent when the time-odd fields of the energy density functional are put 
to zero. Thus, efa is adequate in many practical applications that do not 
require high accuracy. 
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Although for the functionals restricted to time-even fields, the time- 
reversed quasiparticle states a and a are exactly degenerate, this does not 
hold in the general case. Here, the blocking prescription may depend on 
which linear combination of those states is used to calculate the blocked 
density matrix. This point can be illuminated by introducing the notion 
of an alispin^^ which describes the arbitrary unitary mixing of \a and 
Xa- x' a — a Xa + b\a ( | a. | 2 + \b\ 2 = 1). As usual, the group of such uni- 
tary mixings in a 2 x 2 space can be understood as rotations of abstract 
spinors, which we here call alirotations of alispinors. If the time-reversal 
symmetry is conserved, the blocked density matrix becomes independent of 
the mixing coefficients (a, b), that is, it is an aliscalar. In the general case 
where time-reversal symmetry is not dynamically conserved, however, the 
blocked density matrix is not aliscalar. Here, the blocked density matrix 
may depend on the choice of the self-consistent symmetries and the energy 
of the system may change with alirotation. 

The key point in this discussion is the realization that blocking must 
depend on the orientation of the alignment vector with respect to the prin- 
cipal axes of the mass distribution. To determine the lowest energy for 
each quasiparticle excitation, self-consistent calculations should be carried 
out by varying the orientation of the alignment vector with respect to the 
principal axes of the system.^ 01 * 102 ! While in many practical applications 
one chooses a fixed direction of alignment dictated by practical considera- 
tions, it is important to emphasize that it is only by allowing the alignment 
vector to point out in an arbitrary direction that the result of blocked cal- 
culations would not depend on the choice of the basis used to describe the 
odd nucleus. Illuminating examples presented in RefP^ demonstrate that 
the choice of the alignment orientation does impact predicted time-odd 
polarization energies. 

Examples of self-consistent HFB calculations of one-quasiparticle states 
can be found in Refs F" ' 1031 ™ (full blocking) and Refs.ES^Hno] ( EFA ) 

9. Summary and conclusions 

The superfiuid DFT based on self-consistent HFB has already become the 
standard tool to describe pairing correlations in atomic nuclei. Such frame- 
work has been implemented in numerous approaches aiming at a consistent 
description of particle-hole and particle-particle channels, and it is gradu- 
ally replacing a much simpler original BCS theory. This is so, because in 
finite systems like nuclei, spatial dependence of particle and pairing fields 
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has to be properly described, especially in the nuclear periphery of weakly 
bound isotopes. In this respect, the BCS theory and its different flavors 
are manifestly deficient EMU] 

In this study, we aimed at presenting some basics of the local superfluid 
DFT along with several aspects of it related to advanced current appli- 
cations. There are, of course, numerous aspects of the HFB theory that 
we could not cover in this limited overview. First, there have been many 
applications of the HFB theory using finite-range interactions, which im- 
ply nonlocal pairing fields. While they are significantly more difficult to 
treat, they do not lead to ultraviolet divergencies. Based on the current 
description of the limited set of nuclear observables related to pairing, it 
is difficult to judge whether the finite range is essential. In fact, one can 
understand finite-range interactions in terms of regularized local function- 
als. Second, we did not discuss various issues related to the restoration of 
particle-number symmetry. Effects of particle-number nonconservation are 
probably little significant in heavy nuclei, but they may become crucial for 
some observables and in specific systems, like, for example, nuclei with only 
few particles in valence shells. Third, we could not cover subjects related to 
the treatment of pairing in high-spin states where the broken time-reversal 
symmetry precludes the use of the BCS theory. Fourth, pairing correla- 
tions impact nuclear dynamics in a profound way. Recently, there have 
been many exciting developments related to the treatment of small- and 
large-amplitude collective motion in weakly-bound superfluid nuclei. Fi- 
nally, we did not discuss details of the HFB theory applied to the isoscalar 
pairing. This channel becomes essential in nuclei with almost equal num- 
bers of protons and neutrons, and numerous applications of the HFB theory 
to this case exist in the literature, see Refs! 25 ! 26 ! Many of these topics are 
discussed in other contributions contained in this Volume P 
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